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Abstract. Stochastic analysis was applied to observa¬ 
tions of spontaneous behavior in the carnivorous mol¬ 
lusc Melibe leonina. Six behaviors were defined that 
could be easily recognized on inspection and it was 
found that transitions between each of these behaviors 
could be fully described by a first-order random process 
without memory of past behavioral choices. The behav¬ 
iors are organized by frequency of transition into two 
modes, a feeding mode and a resting mode. Transitions 
within modes are more likely than transitions between 
modes, and the feeding and resting modes arc linked by 
a preferred pair of behavioral transitions. The amount of 
time spent in the feeding mode is positively correlated 
with body size, but the average length of a feeding epi¬ 
sode is independent of size. This suggests that body size 
regulates the probability of entry into feeding behavior 
but does not influence the basic pattern of feeding. In the 
presence of food the animals express nearly continuous 
feeding behavior, suggesting that food reduces the prob¬ 
ability of exiting the feeding mode. This model of spon¬ 
taneous behavior in Melibe is used to form hypotheses 
amenable to further exploration through neurophysio¬ 
logical experiments. 

Introduction 

Analysis of animal behavior provides the background 
necessary for studies of the underlying neural circuitry 
and circuit function. Such analysis sometimes permits 

Received 2 December 1996: accepied 28 February 1997. 

PresenI address: Depl. of Zoology. University of Washington, Seat¬ 
tle. W \ 

‘ F It address: Depl. of Neurobiology, Bryan Research Bldg., 
Duke ( ersilN Medical Center, Durham, NC 27710. 

^ Autl ! horn correspondence should be addressed. 


useful predictions to be made about specific properties of 
the underlying neural networks—for example, the extent 
to which stochastic processes in the nervous system in¬ 
fluence the ordering of behaviors in a sequence (Heili- 
genberg, 1973). Molluscs have been used to advantage 
by neuroethologists because the relatively limited behav¬ 
ioral repertoires and accessible nervous systems of these 
animals permit investigation of synaptic connections be¬ 
tween identified neurons (Davis et ai, 1974; Kandel, 
1976; Gelpcrin, 1983; Willows, 1985; Getting, 1989). 
Most studies have concentrated on reflex or conditioned 
responses evoked by stimuli presented in defined pat¬ 
terns, an approach that simplifies the design of electro- 
physiological experiments. Important generalizations 
about the organization of small neural networks have 
emerged from such studies, including the idea that ani¬ 
mal behaviors are organized in hierarchies of succes¬ 
sively more dominant or more strongly commanded re¬ 
sponses (Tinbergen, 1951; Davis et ai, 1977; Getting, 
1989). 

Concentration on evoked and conditioned behavior 
ignores, however, the phenomenon of unstimulated be¬ 
havioral choice (Lorenz, 1981). We recorded the tempo¬ 
ral sequence of behaviors in the nudihranch Melibe Ico- 
nina as it evolved spontaneously in the absence of overt 
stimulation, then used standard techniques of random- 
process analysis to model this behavioral sequence 
(Chatfield and Lemon, 1970). Our results lead to several 
conclusions. (1) In the absence of overt stimuli, the 
choice of behavior is well described by a first-order ran¬ 
dom process. (2) Consideration of the transition proba¬ 
bilities of this process indicates that spontaneous behav¬ 
iors in Melibe fall into two behavioral modes. (3) There 
is a preferred transition that links the behavioral modes. 
(4) The presence of strong stimuli such as food organisms 
influences the individual transition probabilities. Fi- 


418 


BEHAVIORAL MODES IN MFAJBE LEOSISA 


419 


nally, (5) physiological parameters, such as body size, 
also influence behavioral transition probabilities. 

MeUbe is well suited to this study. We were able to 
define a set of six canonical behaviors that arc easily rec¬ 
ognized on inspection. Transitions between these behav¬ 
iors were infrequent and could be adequately scored by 
making observations at 15-min intervals. This permitted 
us to examine behavioral transitions over an extended 
time period and to perform replicate sets of observations. 
Furthermore, the Mclibe central nervous system is suit¬ 
able for microelectrode and optical recording methods 
for monitoring cellular activity in identified neurons in 
restrained whole animal preparations (Willows, 1973; 
Cohen et al., 1991). This provides the opportunity to 
study mechanisms of spontaneous behavioral expression 
at the level of neural networks. 

Materials and Methods 

Melibe teonina, a carnivorous nudibranch of the fam¬ 
ily Tethyidae, is distinguished by an extended oral veil 
used in prey capture (Agersborg, 1923). Two groups of 
10 animals were studied between 1 and 29 May. Adult 
specimens weighing 40-250 g were collected from the 
Maeweystis kelp forest in Monterey Bay, California. The 
animals were held in 45-liter tanks supplied with contin¬ 
uously flowing sand-filtered seawater for 2-28 days be¬ 
fore observations were begun. Changes in motivation 
were minimized by allowing the animals to acclimate to 
the experimental tanks and by simulating a natural envi¬ 
ronment. No correlation was found between the number 
of days in captivity and the behavioral pattern. For ob¬ 
servation, animals were placed in five identical 17-liter 
acrylic plastic tanks provided with flowing sand-filtered 
seawater and a frond of kelp. The tanks rested on sea 
tables in a roofed aquarium building where they were 
protected from direct sun and rain but exposed to natu¬ 
ral cycles of light and temperature. Methylene blue dye 
was injected into the first right ceratum of one animal of 
the two animals in each tank to facilitate identification. 
The animals were fed every other day at the same time 
with the same amount of their natural food, live mysid 
shrimp collected from the wild. The animals were other¬ 
wise left undisturbed for eight days of observation, which 
included four nonfeeding days from each group. 

We define six behaviors that occur frequently but 
without overlap in the absence of overt stimulation (Fig. 
1). Together they account for greater than 90% of the 
activity in the observation tanks. The behaviors are dis¬ 
tinguished by characteristic body postures involving the 
oral hood, the rhinophore processes, and the cerata, as 
well as by the presence or absence of pedal locomotion. 
The six behaviors are as follows: (1) Feeding (F): rhyth¬ 
mic extension and contraction of the oral hood in a cast¬ 


ing motion through the water (Hurst 1968; Watson and 
Trimarchi, 1992). (2) Open hood (OH): full extension of 
the oral hood with the cerata and rhinophore processes 
extended. (3) “Alert'Vprocesscs extended posture (AL): 
oral hood partially closed with the oral tentacles tucked 
inside, cerata and rhinophore processes extended. (4) 
Roaming/open hood (R/OH): OH behavior combined 
with pedal locomotion. (5) Roaming/processes extended 
(R/AL): AL behavior combined with pedal locomotion. 
(6) Resting (RST): locomotion absent, hood closed and 
held against the substrate, cerata contracted against the 
sides of the body, rhinophore processes folded down. 

Several other behaviors occurred too infrequently for 
analysis. These included (1) pedal locomotion in the 
RST posture; (2) the ''crumple’' reflex, an alarm re¬ 
sponse elicited by tactile or vibratory stimuli and charac¬ 
terized by forceful contraction of most of the body wall 
musculature (Scott, 1990); (3) swimming; (4) egg-laying. 
Copulation was observed occasionally but was not 
scored independently because the six behaviors we fol¬ 
lowed occurred at the same frequency independent of 
copulation. 

Observations were made at 15-min intervals for 6 
hours a day beginning between 0800 and 1200 h. During 
preliminary studies, we observed behavior continuously 
for periods of several hours. By parsing the data into time 
periods of varying duration, we determined that a 15- 
min sample interval is adequate to capture behavioral 
transitions. Observations from one nonfeeding day for 
the 10 animals in group I are tabulated in Figure 2. It 
is apparent that behavior changes infrequently and that 
individuals tend to perform the same behavior over sev¬ 
eral time points. On feeding days the animals spent 
about 80% of their time performing feeding behavior (F), 
indicating that the presence of prey organisms has a pow¬ 
erful influence on behavioral choice (Watson and Tri¬ 
marchi, 1992; Watson and Chester, 1993). Analysis was 
limited to nonfeeding days in order to include a more 
diverse range of behaviors and to obtain observations un¬ 
der relatively constant physiological conditions. 

Results 

The frequency of occurrence of each of the 36 possible 
behavioral transitions, including consecutive occur¬ 
rences of the same behavior, was tallied by observing be¬ 
havior at 15-min intervals. The data are expressed in Ta¬ 
ble 1 as transition probabilities for each of two replicate, 
4-day experiments. The sum of the transition probabili¬ 
ties from a given starting behavior is normalized to 1.0. 
The predominance of certain one-step transition proba¬ 
bilities over others suggests that the six behaviors can be 
grouped into two distinct behavior modes (Fig. 3). Be¬ 
haviors F, OH, and R/OH form the feeding mode, while 
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Figure I. Body postures during six behaviors: (A) Open hood (OH) and roaming/open hood (R/OH) 
behaviors. (B) Resting (RST) behavior. (C) Feeding (F) behavior. (D) “Alert’Vprocesscs extended (AL)and 
roaming/alert (R/AL) behaviors. 


behaviors RST, AL, and R/AL form the resting mode. 
This grouping of behaviors is based on the observation 
that transitions among mode elements occur much more 
frequently than transitions to behaviors outside of the 
mode. The frequencies of intermode transitions are 
listed in Table 11. This qualitative observation is borne 
out by calculating between-mode transition frequencies 
for all possible ways of grouping the six behaviors into 
two modes. The correct mode structure will give fewer 
between-mode transitions than any other grouping. 

There are 25 ways to structure two modes: 10 in which 
each mode contains three behaviors, and 15 in which one 


mode contains two behaviors and the other contains 
four. For data set 1, the number of between-mode tran¬ 
sitions for the mode structure proposed in Figure 3 was 
145, whereas between-mode transitions for the other 
possible groupings ranged from 267 to 389. This sup¬ 
ports the behavioral relationships diagrammed in Figure 
3. According to this analysis, animals stay in one mode 
for an average of 1 hour and 11 minutes (5.2 observation 
intervals) before switching to the other. 

Transitions between modes occur most frequently be¬ 
tween behaviors OH and AL in either direction. The OH 
to AL transition accounts for 35% of all transitions from 
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Figure 2. Sample data sei from one full day ot observations. Observations of each of 10 animals were 
made at l5-min intervals for 6 h. The duration of a behavior in continuous intervals is referred to as its 
dwell time. An example of a seven-interval dwell lime in the open hood (OH) behavior can be seen in 
animal 4A from intervals six to twelve. Abbreviations: Feeding (F), open hood (OH), “Alert’Vprocesses 
extended posture (AL), roaming/open hood (R/OH), roaming/processes extended (R/AL), resting (RST). 
crumple (CR), and swimming (S). 
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Table I Tabic 11 

Transition prohahiliucs between all behaviors for both sets of data Transition probabilities within and between modes 


SET I 

Ending behavior 

Slanting 


behavior 

F 

OH 

R/OH 

AL 

R/AL 

RST 

F 

0.597 

0.223 

0.095 

0.038 

0.019 

0.028 

OH 

0.183 

0.486 

0.115 

O.lOl 

0.014 

0.101 

R/OH 

0.299 

0.253 

0.218 

0.080 

0.151 

0.226 

AL 

0.041 

0.142 

0.027 

0.473 

0.122 

0.196 

R/AL 

0.085 

0.149 

0.085 

0.298 

0.149 

0.234 

RST 

0.040 

0.081 

0.054 

0.117 

0.054 

0.664 

SET 2 







Starting 

behavior 



Ending behavior 



F 

OH 

R/OH 

AL 

R/AL 

RST 

F 

0.459 

0.208 

0.138 

0.101 

0.075 

0.019 

OH 

0.188 

0.370 

0.087 

0.268 

0.043 

0.043 

R/OH 

0.293 

0.200 

0.227 

0.093 

0.120 

0.067 

AL 

0.065 

0.094 

0.039 

0.526 

0.101 

0.175 

R/AL 

0.135 

0.122 

0.i35 

0.446 

0.122 

0.041 

RST 

0.005 

0.020 

0.010 

0.256 

0.054 

0.655 


Each row of probabilities totals 1.00. Probabilities for transitions oc¬ 
curring within a behavioral mode are given in boldface. 


the feeding to the resting mode, and the AL to OH tran¬ 
sition accounts for 30% of transitions from the resting to 
the feeding mode. None of the 16 other possible bc- 
tween-mode transitions accounted for more than 15% of 
the transitions in either direction. 

/ndividual variahilily 

The time spent in each behavior varied a great deal 
between individual animals. As a measure of this vari- 


FEEDINGMODE ^ 17 ^ 

OH R/OH 

-----it.-■ 

AL < > R/AL 

RESTING MODE RST 

Figure 3. Eeeding and resting modes. Arrows indicate the 10 most 
frequent iransitions between different stales. The mode structure min¬ 
imizes intermode transitions: all other possible mode arrangements 
were considered (see lext). 


Probability that 
the next transition is to: 


Feeding mode Resting mode 

Starting from: 

Feeding mode 0.79 0.21 

Resting mode 0.18 0.82 


ability, wc calculated coelficients of variation (CV) for 
the number of transitions from a behavior back to itself 
{i.c., F F) using the second set of 10 animals. The CV 
ranged from 64% to 157%. Some of the variability can be 
linked to a physiological correlate. For example, we 
found that body weight and time spent in the feeding 
mode are positively correlated (/• = 0.87: Fig. 4). An 80-g 
individual spent 10% of the time in the feeding mode as 
compared to 80% for a 250-g animal. A similar correla¬ 
tion exists between size and time spent in the feeding (F) 
behavior within the feeding mode. The 250-g animal 
spent 58% of the total time in F, but the 80-g animal did 
not exhibit this behavior at all. This bias on both inter- 
and intramodal transition probabilities may reflect a 
difference in the nutritional needs of small vs. large ani¬ 
mals and is worthy of further study. We did not find a 
correlation between body size and average dwell time in 
the F behavior in the absence of food organisms. This 
suggests that large animals enter feeding behavior more 
often than small ones, but that the duration of feeding 
episodes is about the same once feeding begins. Although 
body size has little effect on feeding dwell time, this pa- 



Figure 4. Body size and total time spent in the feeding mode are 
correlated. Data taken from set 2.1 olal time in the feeding mode from 
Figure 3 was calculated by counting the total number of intervals spent 
in feeding mode behaviors. Body volume was determined by w’ater dis¬ 
placement. The correlation coefficient is r = 0.87. 
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rameter is influenced > i!.e presence of food; Watson 
and Trimarchi (1'^' ' c shown that dwell time varies 

inversely with the , entration of available prey. 

Two-siep trc’! piobabililies 

A sequenre ol behavioral acts is stochastic if the tran¬ 
sition from one behavior to the next is independent of 
preceding behavioral changes. If this is so, the frequency 
of two-step transitions observed experimentally should 
match the transition frequencies predicted from a two- 
step stochastic model. Expected probabilities for two- 
step transitions were calculated from the measured one- 
step probabilities following the equation /"(a -* b ^ c) 

= F(a b) * P{b -* c). where a. b. and c represent any 
of the six canonical behaviors. Behavioral observations 
from the 10 individuals in each data set were pooled and 
used to tabulate the probability of occurrence of each of 
the 216 possible two-step behavioral transitions. Transi¬ 
tions from four of the six possible starting behaviors oc¬ 
curred frequently enough to permit this extensive analy¬ 
sis (F. OH, AL, and RST). The expected and observed 
probabilities were placed in 6 X 6 matrices for compari¬ 
son. The matrices for starting behavior a = AL from data 
set 2 are tabulated and plotted in Figure 5. The expected 
and observed three-dimensional plots of behavioral fre¬ 
quency show peaks and valleys at the same positions and 
with approximately the same amplitudes. 

We attempted to test the significance of differences be¬ 
tween the expected and observed 6x6 matrices by cal¬ 
culating chi-square values. When the two data sets were 
treated separately, chi-square scores ranged from 42.5. 
for set 2, starting behavior AL, to 62.9 for set 1, starting 
behavior F. When both sets of data were pooled, chi- 
square scores ranged from 66.5 to 81.5. All the chi- 
square values were significant by standard tests using 25 
degrees of freedom (/’ < 0.001 Ko P < 0.01). The chi- 
square statistic is known to be unreliable, however, in the 
case of behavior-sequence matrices with cells containing 
low values (Chatfield and Lemon, 1970), and the two- 
step transition matrices we tabulated have a number of 
cells with values ofO. Also, variaii sn in transition rates 
was high among individuals (CV lor F — F as high as 
157%), which adds to the chi-square score. The presence 
of this type of added uncertainty is demonstrated by the 
fact that higher scores result when data from two periods 
of observation are pooled. Instead of a chi-square test, 
Chatfield and Lemon (1970) suggest selecting the most 
frequent sequences and comparing expected and ob¬ 
served values by inspection. The plots of Figure 5 follow 
that principle and show that the two-step transition pre¬ 
dictions match the observations fairly well. 

Behavioral dwell-tbne distributions 

We defined behavioral dwell time as the duration of 
a continuous single bout of the same behavior. This is 


equivalent to the 'interval distribution’ (Heiligenberg, 

1973) or the 'hazard function’ (Cox and Lewis, 1966) but 
differs in that exits into any other behavioral state are 
allowed. Dwell-time histograms for four of the six behav¬ 
iors are plotted in Figure 6. If behavioral choice is deter¬ 
mined stochastically, the distribution of dwell times will 
be exponential. R/OH and R/AL behaviors did not oc¬ 
cur frequently enough to permit dwell-time analysis. OH 
and AL dwell times were exponentially distributed. F 
dwell times were also exponentially distributed, with the 
exception of a few lengthy episodes. The presence of 
longer episodes has implications for the organization of 
the feeding neural network and suggests that the network 
may be modulated by stimuli such as the presence of 
stray food organisms or by a long-lasting effect of food 
presented on the previous day (Watson and Chester, 
1993). The RST dwell-time histogram was best fit by the 
sum of two exponentials. This suggests two possibilities: 
first, there may be two resting states which have indepen¬ 
dent transition probabilities but are indistinguishable 
by observation; alternatively, the transition probabilities 
for exiting RST may vary with time spent in the RST 
behavior. 

Discussion 

These observations demonstrate that spontaneously 
evolving sequences of behavioral transitions in Melibe 
leonina are stochastic such that the probability of a given 
transition is independent of preceding ones. At the same 
time, behavior is organized by the fact that transitions 
between behaviors within a mode are more likely than 
transitions between behaviors in different modes. This 
creates a cul-de-sac w ithin the transition-probability ma¬ 
trix so that one behavior pattern can only be reached by 
way of another. As a result, the organization is not in¬ 
consistent with hierarchical models of behavioral orga¬ 
nization like that proposed by Tinbergen (1951) since 
specific stimuli might move the animal into a mode in 
which certain behaviors are more likely than others even 
though the individual transition probabilities are con¬ 
stant and independent. This view is supported by evi¬ 
dence that the neural circuits underlying spontaneous 
behavioral choice are influenced by external cues such as 
the presence of prey and by physiological signals corre¬ 
lated with body size. 

Spontaneous behavior in Melibe demonstrates two 
features characteristic of a stochastic system. First, the 
two-step transition probabilities predicted from consid¬ 
eration of measured single-step transition frequencies 
closely match the two-step frequencies observed experi¬ 
mentally. This shows that behavioral transitions are in¬ 
fluenced by the starting behavior but not by previous be¬ 
haviors. Second, for three of the four behaviors for which 
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Figure 5. Top. Expected and observed Iwo-siep transitions with data from starting behavior AL (set 2). 
Expected numbers were calculated from single-step transition frequencies using the equation P(a b -► 
c) = /’(a)/^(a bV^(b -► c) (see text). Bottom. 3-D plots corresponding to the expected and observed two- 
step transition matrices. 


there was sufficient data to evaluate behavioral dwell- 
time distributions, dwell times were distributed expo¬ 
nentially. This supports a model based on independent 
and unvarying transition probabilities. The dwell-time 
distribution for the exception, the RST behavior, was 
best fit by the sum of two exponentials. This result might 
be explained by the existence of two resting states with 
different average duration. We were unable to identify 
multiple resting states but, like sleep and quiescent be¬ 
havior in other animals, the two states may have a 
different neurophysiological basis while appearing iden¬ 
tical to the observer. Another possible interpretation is 
that the exit probability from RST decreases with time 
spent in the RST state, perhaps reflecting a decrease in 
arousal. The data are not sufficient to evaluate these pos¬ 
sibilities, but they raise questions of a general nature that 
might be accessible to neurophysiological investigation 
in Melihe. 

A properly that emerges from a consideration of the 
transition-probability matrix is that spontaneous behav¬ 
ior can be interpreted as a series of jumps between two 


modes named ihe feediiig mode and the resting mode. 
This grouping of behaviors follows directly from the one- 
step transition probabilities and the observation that 
transits within modes are far more probable than transits 
between modes. It is a robust conclusion based on over 
1900 observations of transitions, far more than the 360 
(10 m) required to analyze a repertoire of n behaviors 
(Fagen and Young, 1978). This organization represents 
a way in which higher order structure can be generated 
by a coupled system of first-order random processes. 

The probability of transitions between modes did not 
fluctuate during 6 h observations, but this period is short 
relative to the life span of the animal. Preliminarx' data 
suggest that transition frequencies do change with matu¬ 
ration and that juveniles demonstrate both a dift'erent 
repertoire of behaviors and different frequencies of tran¬ 
sitions between the behaviors they share with adults (S.E. 
Gelber, pers. comm.). In addition, Watson and Trimar- 
chi (1992) have shown that the presence of prey organ¬ 
isms, sensed in response to chemical and tactile cues 
(Watson and Chester, 1993), intlucnces feeding dwell 








424 


A. E. SCHIVELL ET AL 



o> 200 r 

1 ■" i 1 I 

150 ^ 

4> 

RST 

L. 

” 100 ■ 

\ 

£ 

UK 


S 50 ■ 

o 

0 ■ 

n 

w , 

1 


0 60 120 180 240 300 

time (minutes) 




Figure 6. Dvvell-time histograms for feeding (F), open hood (OH), resting (RST), and alert-processes 
extended (AL) behaviors. The OH and AL distributions are fit by exponentials (solid curve), consistent with 
the exit from these behaviors being random. The feeding (F) dwell-time distribution is also approximately 
exponential, except for outlying points that account for l%of the data. The outlying points with long dwell 
times may reflect behavior driven by external stimulation. Two fits are plotted for the resting (RST) data, a 
single exponential (solid line) and the sum of two exponentials (dashed line). All distributions were fitted 
using a nonlinear least-squares method. 


time in a fashion that is graded with the concentration of 
prey, 

A particular pair of behavioral transitions, OH^AL, 
acts as a gateway between the feeding and resting behav¬ 
ior modes. The OH and AL behaviors share features of 
body posture such as the positions of the rhinophores 
and cerata, but they different in oral hood posture. The 
change between behavior modes, therefore, is bridged by 
a relatively small postural adjustment, but it signals a 
long-term change in the animal's behavior pattern. A 
neurophysiological interpretation is that slight but mea¬ 
surable shifts in the activities of neurons involved in OH 
and AL behaviors can result in a stable, lasting shift in 
motor-system function. 

Increased body size is correlated with an increase in 


time spent feeding, but the average dwell time of feeding 
episodes is independent of size. This suggests that the 
probability of entering feeding behavior increases with 
body size. Conversely, it is apparent that the presence 
of food organisms decreases the probability of exit from 
feeding but does not have a clear effect on entry into feed¬ 
ing behavior. If we assume that separate entry and exit 
circuits govern transitions into and out of the feeding 
state, it should be possible to measure changes in excit¬ 
ability in those circuits that correspond to the behavioral 
observations. One simple and testable prediction is that 
body size modulates the entry circuit, while the presence 
of food modulates the exit circuit. Modulation of the exit 
circuit by food could confer competitive advantage by 
maximizing foraging efficiency in the face of a variable 
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food supply. Such shifts in behavioral transition proba¬ 
bilities might be mediated by neuromodulatory or hor¬ 
monal effects that cause a dynamic restructuring of neu¬ 
ral networks by influencing either postsynaptic excitabil¬ 
ity or the amplitudes of synaptic potentials (Getting, 
1989; Katz 1994). 

Neuronal circuits that mediate behavior in MeliheQ^LVi 
be studied with intracellular recording in a whole-animal 
preparation using techniques pioneered by Willows 
(1973) and Getting (1981), Using this preparation, we 
have made progress in elucidating the networks involved 
in swimming locomotion and feeding movements of the 
oral hood; Cohen et ai (1991) have succeeded in using 
optical imaging techniques to monitor activity simulta¬ 
neously in dozens of neurons in the Afelibe buccal gan¬ 
glion. Neurophysiological methods are available, there¬ 
fore, to study changes in neuronal function that are cor¬ 
related with behavioral choice. The behavioral model 
developed here provides a framework for the interpreta¬ 
tion of physiological experiments. 
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